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The sufficient conditions for existence and uniqueness of continuous solutions of the Volterra 
operator equations of the first kind with piecewise continuous kernel are derived. The asymptotic 
approximation of the parametric family of solutions are constructed in case of non-unique 
solution. The algorithm for the solution's improvement is proposed using the successive approximations 
^ method. 
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1. Introduction 

Let us introduce area D = {s, t; < s < t < T} in the plane s,t and define continuous 



00 

Q 

functions s = ojj(t),z = l,n, which have continuous derivatives for t G (0, T). Let us suppose 
that oji(O) = 0, < ai(t) < ■■■ < a n -i(*) < t for t G (0,T), < a[(0) < ■■■ < a'^O) < 1, 
and functions s = a,(i), % = 0,n, cto(t) = 0, a n (t) = t, split the area D into disjoint areas 

n 

D 1 = {s,t : < s < cti(t)}, Di = {s,t : oti_i(t) < s < aii(t),i = 2,n}, D — (J A- Let 

O 
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us introduce biparametric family of linear continuous operator-functions Ki(t,s), defined for 



t, s G Di, i = 1, n, which are differentiable wrt £ and acting from Banach space £?i into Banach 
space E 2 . Therefore Ki{t,s) G C{E X -> E 2 ) and G -> S 2 ) for t,s e~D~,i = l~n. 

Let the space of continuous functions x(t) defined on [0,T] with ranges on E\ be denoted as 
C([o,t];£i)- Let us introduce the integral operator 

t n Oi(t) 

I K{t,s)u{s)ds = y y [ Ki(t,s)x{s)ds (1) 



i=1 ai-i(t) 



with piecewise kernel 



K^s), t,se Di, 

K(t,s) = { (2) 

K n (t,s), t,s G £>„• 
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Academic Staff of Innovative Russia" within project No. 2012-1.2.2-12-000-1001-012. 
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In this text we concentrate on the equation 

t 

J K(t,s)x(s)ds = f(t), (3) 
o 

where function f(t) and function f'(t) with ranges in E 2 are defined and continuous for t G 
[0, T], /(0) = 0. Let us call the equation (3) as the Volterra operator equation of the first kind. 
The objective is solution to the equation (3) in C([o,t];£i)- Such problem and its numeric solution 
for Ei = E2 = K 1 , n = 1 has been studied in number of papers, see e.g. [1-4]. Differentiation of 
the equation (3) lead to the integral-functional equation and its solution is not unique in the 
general case. The reader may refer to [5]. Therefore construction of the solution to equation 
(3) can not be carried out with classic analytical mathods of the Volterra equations [6]. In this 
paper we address the equation (3) in the general case using the theory of operator equations 
with functionally perturbated argument of neutral type [7, 8]. It is to be noted that the theory 
of the Volterra operator equations was initiated in [9]. The principal results in this area with 
applications in integral geometry have been obtained in scientific school of M.M.Lavrentiev (see 
there references e.g. in the monograph [10]). 

The paper is organized as follows. In Section 2 we continue our studies on the Volterra 
equations presented in papers [5, 11-17] and provide the sufficient conditions of existence and 
uniquness of continuous solution of the equation (3) with piece-wise continuous kernel (2). The 
special case of such equations are the Volterra equations which model the evolving dynamic 
systems [1-3]. To the best of our knowledge such equations have not been studied in the 
literature. In the Section 2 the desired solution is constructed based on the "step method" [8] 
from the theory of functional equations combined with the method of successive approximations. 
In Section 3 and Section 4 we address the most interesting case when equation (3) has the family 
of solutions which depends on the free parameters. We propose the method for construction of 
asymptotic approximations of parametric solutions and we design the method of the solution's 
refinement using the successive approximations. 

In this work it is supposed that an operator K n (t,t) has limited inverse operator for 
t G [0, T]. Norms \\Ki{t, s)\\c(e 1 ~>e 2 )-, 1 1 9Ki ^' s ^ \ \c(e 1 ^-e 2 ), * — l; n °f the linear operators are 
defined and they are continuous functions for t, s G Di, t G [0,T]. The conventional notations 
and terminology from [19, 20] are employed. 

2. Sufficient Conditions for Existence of the Unique Continuous Solution 

Since /(0) = differentiation of equation (3) leads to the equivalent functional-operator 
equation 

n— 1 s 

F(x) ^ K n (t,t)x(t) +Y t a'i(t){Ki(t,a i (t))- 

i=i ^ 
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-K i+1 (t, Oi(t)) \x(oi(t)) I K^(t, s)x(s) ds - f'(t) = 0, 

J i=l /,n 



(4) 



£*i-i(t) 



where a = 0, = t. 

Let us introduce the function 



Tl-l 

def 



i=i 

where || • ||£(_Ei^Si) is the linear operator's norm. Let the following condition be fulfilled: 
(A)L>(0)<1, sup \\K-\t,t)K{t,s)\\ c{ E 1 ^E 1 )<c<oo. 

0<s<t<T 

Clearly the inequality D(0) < 1 is fulfilled if |a^(0)| are sufficiently small. Here and below 

n 

an operator K(t, s) is defined with formula (2) in U-Dj. It's derivative wrt t in common sense 

l 

is defined as 

f KPfas), t,seD 1 , 

K^\t,s)=l (5) 

[ K^\t,s), t,se D n . 

n 

for t,s G (J A- 
l 

Theorem 1. (Sufficient conditions for existence and uniqueness of the solution) 

Let condition (A) be fulfilled, all the operators Ki(t, s) in (2) are continuous, and wrtt they also 
have continuous derivatives, vector f(t) has continuous derivatives, /(0) = 0. Then equation (3) 
has unique solution in the class of continuous functions C([ 0i t];Ei)- And moreover the solution 
can be constructed using the step method combined with successive approximations. 

Proof. Let us apply the operator K~ l (t,t) to the equation (4). We get the following equation 

x(t) + Ax + Kx = J(t), (6) 

where the following notations are introduced 

A(t)x = K-\t,t) n Y j af\t){K t {t,a l {t)) - K l+1 (t,a t (t)))^t), 

1=1 

Kx = £ jf K- l {t,t)K i l ) {t ) s)x{s)ds, J{t) = K-\t,t)f^(t). Let us fix q < 1 and select 

i=lQi_i(t) 

hi > such as max D(t) = q < 1. Such hi > can be found because of condition (A), and 
0<t<hi 

since the operator-functions Ki(t, s) are continuous in operator topology and because functions 
cti(t) are continue with their derivatives. Lets < h < min{/ii, where constant value c is 

defined in the condition (A). Let us split [0,T] into 

[0, h], [h, h + eh], [h + eh, h + 2eh], .... (7) 
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The contraction of the desired solution x(t) onto the interval [0,h] we denote as x (t), the 
contraction onto 

I m = [(1 + (m- l)e)h, (1 + me)h], m = 1, 2, . . . 
we denote as x m (t). Let us select e from (0, 1] such as for t G I m the "perturbed" arguments 



m— 1 



^ |J i = 1, n — 1. If < a- 1 ^) < yq— for t G [0,T), i — l,n — 1, then the mentioned 



fc=i 



l+e 



above inclusion is fulfilled on [0,T). Such inclusion enable the application of the "step method" 
(readers may refer to [8], p. 199)) from the theory of functional-differential equations for 
construction of the solution x(t). 

For construction of the element x (t) G C^ 0th ] :El ) we construct the sequence {^(t)} : 

xZ(t) = -Axr 1 -Kxt- 1 + f(t), 
x° (t)=J(t),te[o,h]. 

Because of selection of h we have an estimate \ \A + K\\c(c {[0 h] . El) -^c {[0 h] . El) ) < 1- 

Therefore for t G [0, h] exists the unique solution Xo(t) of the equation (6). The sequence 

£q(£) converge uniformly to this solution. Let us continue the process of the desired solution 

construction for t > h, i.e. on 7 n , n — 1, 2, For sake of clarity let e — 1 in (7). 

Then once we have the element x (t) G C([o j h] ; s 1 ) be constructed we search for the element 

Xi(t) in the space < C{[h,ih\;E 1 ) of continuous vector-functions. Let us find Xi(t) from the Volterra 

equation of the 2nd kind 



(t) + J K- 1 (t,t)K> t (t,s)x(s)ds = 



h 

h 



= f(t) -Ax - J K~\t, t)K' t {t, s)x (s) ds 





using the successive approximations. Here x (h) = Xi(h). 
Let us introduce the continuous function 

UM, o<t<ft, 

\ xi(t), h<t<2h, 

which is the contraction of the desired solution x(t) onto [0, 2h]. Then the element X2(t) G 
C([2/ l) 3/ l ];£ 1 ) can be computed with successive approximations from the Volterra integral equation 
of the second kind 



(t) + J K- 1 (t,t)K^s)x(s)ds 



x v 

2h 

2h 



= f(t)-Ax 1 - J K- 1 (t : t)K' t (t : s)x 1 (s)ds. 
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The desired solution x(t) G C([ 0i t] ; e 1 ) to the equation (3) can be constructed on iVth step, 
N>1 



Example 1. 
Integral equation 

t/2 t 

J K ± (t - s)x(s) ds + J K 2 (t - s)x(s) ds = f(t), < t < T, 

t/2 

/(0) = 0, Ki(t — s) — K 2 (t — s) + E, Ki, K 2 are matrices m x m, E is the unit matrix, 

1-^2 1 (0)|£(K m -i>IR m ) < 2, 

matrix K 2 (t) and vector-function f(t) = (fi(t), . . . , f m {t))' have continuous derivatives wrt t, 
and satisfy conditions of the Theorem 1 and has unique continuous solution. 

Example 2. 
Boundary problem 

t/2 / \ t 

I (*$fi + x(t,y)) dt +fj^dt = f(t,y), < t < T, < y < 1, 
x(t,0) = 0, x{t,l)=0, 

where function f(t,y) is continuous wrt y and has continuous derivative wrt t, f(0,y) = 0, 
satisfy the conditions of the Theorem 1. Desired continuous solution can be constructed as 
solution to the equivalent equation 

t l 

1 

xC 

2 

o o 
with contracting integral operator where 

' (f - i)y, y<Z, 
(y-i)C, Z<y 

based on the successive approximations method. 



t i 
'(t,y) = ~ J G(y,04t,0dt + J G(y,t)fl(t,t)d£ 



G(y,0 = 



3. Construction of the asymptotic approximation x(t) of parametric family of solutions 

Let the following condition be fulfilled 

TV 



(B) Exist operator polynomials V\ — K iv ^t v s^, % — where K^ v ^ £ C(Ei — > E 2 ), 

V+fJ,= l 

N 

are linear continuous operators, vector-function f N (t) = fv? ' , polynomials af(t) = 

v=l 



N 

J2a iu t u ,i = — where < an < a 2 i < «2i < ••• < oe n -i,i < 1) are sucn as 
for t — > +0, s — > +0 the following estimates are fulfilled \\Ki(t,s) — Vi(t, s)\\c(e 1 -^e 2 ) = 
0((t+s) N + l ), < = M, \\f(t)-f N (t)\\ E2 = 0(t N+1 ), \a t (t)-af(t)\ = 0(t N+1 ),i = 1^1 

Expansion of powers t, s in the condition (B) we call as "Taylor polynomials" of the corresponding 
elements. 

Let us introduce j-parametric family of linear operators 

n-l 

B(j) = K n (0, 0) + ^K(0)) 1+ ^(0, 0) - K i+1 (0, 0)), 
i=i 

j £ [0, oo). Operator B(j) which corresponds to the main "functional" part of the equation (4). 
We denote operator B(j) as characteristic operator of equation (4). 

Let us follow the section 2 and consider equation (4) which is equivalent to the equation (3). 
In contrast to the section 2 here we do not suppose that homogenius equation for equation (3) 
has only the trivial solution. Therefore solution to the integral-functional equation (4) can be 
non-unique. We follow [5] and search an asymptotic approximation of the particular solution 
of non-homogenius equation (4) as following polynomial 

N 

x{t) = Y. x ^ t3 - ( Q ) 

3=0 

Let us demonstrate that coefficients Xj with ranges in E\ in the general non-regular case 
depends on lnt and free parameters. That is in line with that fact of existence of nontrivial 
solution of the homogeneous equation. 

There are regular and non- regular cases when it comes to determination of the coefficients Xj. 

Definition 1. j* is the regular point of operator B(j), B(j*) has bounded inverse operator 
and non-regular point otherwise. 

3.1 Regular case: characteristic operator B(j) has bounded inverse operator for 

j e(o,i,...,N) 

In this case the coefficients Xj are constant vectors from E\. Indeed, let us substitute (9) 
into the equation (4). Then using the method of successive approximations and condition (B) 
we get the recurrent sequence of linear equations wrt Xj : 

B(0)x = /'(0), (10) 

B(j) Xj = Mj(x , . . .,^-0, j = 1, . . . ,N. (11) 
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Vector Mj can be presented via solutions x , . . . ,Xj-i and coefficients "Taylor polynomials" 
from the condition (B). 

Since operators B(j) are invertible in the regular case then vectors Xo,...,Xn can be 
uniquely defined and therefore the asymptotic (9) can be constructed. 

3.2 Non-regular case: operator B(j) for j e (0,1,..., N) has non-regular points 

Let us introduce definitions: 

Definition 2. Value j* is simple singular Fredholm point of the operator B(j), if B(j*) is 
Fredholm operator ([20] p. 219), det[< B^(j*)<f)i,if) k >Y ik=1 ^ 0, where is basis in 

N(B(j*)), {ipijl is basis in N(B'(j*)), B'(j*) is conjugate operator, B^(j) is derivative wrt j, 
computed with j = j*. 

Definition 3. Let B(j*) be Fredholm operator, j* we call singular Fredholm point of index k + 1 
ifN(B(j*))c f]N(B^(f)), 



i=i 



r 

det <B (k+1 \f)&^ k > 



^0, k > 1. 



i,k=l 



n-1 



It is to be noted that B^(j) = £ (o4(0)) 1+ %f (K^O, 0) - K i+1 (0,0)), where a { = lnaj(0). 

i=i 

Remark 1. According to the definition 3 index of the simple Fredholm point is 1. If E\ = E2 = M 1 
then B(j) is function of j. In this case the definition 2 means that j* is single root of the equation 
B(j) = 0, definition 3 means that j* is (k + l)-multiple root of this equation. 

Let us demonstrate that in non-regular case the coefficients Xj will be polynomials of power 
In t and depends on arbitrary constants. Order of polynomials and number of arbitrary constants 
are connected with indices of singular points of the operators B(j) and dimension N(B(j)). 

Indeed since coefficient x in the nonregular case can depends on lnt, then based on the 
method of undetermined coefficients x can be determined as the solution to the difference 
equation 

n-1 

K n (0, 0)x o (z) + ^(0)(Ki(0, 0) - K i+1 (0, 0))x o (z + a,) = f'(0), (12) 
i=i 

where a, = \na'(0),z = Int. 

Here there are three cases: 
1st Case. 

Operator B(0) has bounded inverse operator. Then coefficient Xq does not depends on z and 
can be determined uniquely from (10). 
2nd Case. 

Let j = be the simple Fredholm point of the operator B(j). Let us search the coefficient xq(z) 
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from difference equation (12) as linear vector-function 

x (z) = x 01 z + x 02 . (13) 
We substitute (13) into (12) and we get the following two equations for determination of vectors 

B(0)x 01 = 0, (14) 
B(0)x 02 + B^\0)x 01 = f'(0). (15) 

r 

Let {4>i}i be basis in N(B(0)). Then x i = Yl Ck4>k- Vector c = (ci, . . . , c r )' can be determined 

k=l 

uniquely from the system of linear algebraic equations 



£ < #«(O)0 fe ,^ > c k =< /'(0), >, i = 1, 



fe=i 



with nonsingular matrix. Further the coefficient X02 can be determined from the equation (15) 
with accuracy up to span(0 1; . . . , <p r ). For this reason we employ the formula 

r 

x 02 = Y,dk<Pk + W(0)-B {1) x 01 ), 
k=i 

where 

r 

T = (B(0) + Y,<;lk> zk)- 1 

k=l 

is Trenogin's regulariser (see [20], p. 221), di,...,d r are arbitrary constants. Therefore in the 
2nd case the coefficient x (z) is linear wrt z and depends on r arbitrary constants. 
3rd Case. Let j = be the singular Fredholm point of the operator B(j) of index k + 1, k > 1. 
The solution x (z) to the difference equation (12) we search as polynomial 

x (z) = x 01 z k+1 + x 02 z k H h x ok+1 z + x 0k+2 . (16) 

Let us substitute of the polynomial (16) into the system (12) and let us take into account the 
following identity 

^-B{j) = £(^(0)) 1+ ^(0,0) -K i+1 (0,0)), 
aj i=i 

where = lno;-(0). Let us match the coefficients on powers of z k+1 ,z k , . . .,z,z° to zero. As 
result we get the sequence of linear operator equations wrt x 0i , x 02 , ■ ■ ■ , £ofc+2 : 
B(0)x 01 = 0, 



£(0)x 2 + 5 (1) (0)( k+ k l )x O i = 0, 
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B(0)x 0k+2 + B^ k+1 \0)x 01 + BW(0)x 02 + ... B^(0)x ok+1 = /'(O). 
Let us follow the definition 3 here. Therefore we have the following inclusion 

W( B (0))cf>(^ 



(17) 



1=1 



J=0- 



Hence 5^(0)%+! = 0, i = 0, k and coefficients x ±, . . . , :rofc+i can be determined from homogeneous 



equation B(0)x = based on formulas x 0i = Cij4>ji i — l,k + 1. Therefore the equation (17) 
becomes 



B(0)x , k+2 + B( k+1 \0)x 01 = f'(0). 



Since B(0) is the Fredholm operator and det 

.1 4£. f 



< (0)&,^ fc > 



(18) 

7^ 0, then vector 



i,k=l,r 



c' = (cn, . . . , Ci r )' can be determined uniquely from the conditions of resolvability of the 
equation (18). Therefore we have 



X0k+2 — Cfe + 2j0j + Xk+2, 

3=1 



Xk+2 is particular solution to the equation (18). Vector c k+2 = f (ck+2,1, ■ ■ ■ ,c k +2, r )', as we U as 
c 1 = (cji, . . . , Cir)' , i = 2,k + 1, remains arbitrary. As result, the 3rd Case the coefficient Xq(z) 
is polynomial of k + 1th power wrt z and depends on r x (k + 1) arbitrary constants. 

Let us apply the method of undetermined coefficients and take into account the equality 

/ tnn k tdt = v +1 Y\-iy k{k - 1 \ V ' ' { \ ~ (g ~ 1)} in fc -t. 



s=0 



We can construct the difference equations for determination of the coefficient X\(z) (z = lnt) 
and next coefficients of the asymptotic approximation (9). Indeed let us take into account the 
definition of the operator F (see equation (4)). Then we get 



F(x) 



x=xo(z)+xi(z)t 



n-1 



K n (0, 0) Xl (z) + ^K(0)) 2 (^(0, 0)- 



(19) 



i=i 



-K i+1 (0,0))x 1 (z + a i ) + P 1 (x (z)) 



t + r(t), 



with estimate r(t) = o(t). Here P\(xq(z)) is certain polynomial wrt z. Its power is equal to 
index of singular Fredholm point j — of the operator B(j). From (19) and because of the 
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estimate r(t) = o(t) for t — > we conclude that the coefficient Xi(z) have to satisfy the following 
difference equation 

n— 1 

K n (0, 0) Xl (z) + ^(a'(0)) 2 (^(0, 0) - K i+1 (0, 0)) Xl (z + <n)+ (20) 
i=i 

+Pi(x (^)) = 0. 

If j = 1 is regular point of the operator B(j) then equation (20) has the solution X\(z) as 
polynomial of the same order as index of the singular Fredholm point j — of the operator 
B(0). If j = 1 also singular Fredholm point of the operator B(j) the the solution can be 

constructed as polynomial of power ko + ki, where ko and k\ are indices of the corresponding 
singular Fredholm points j = 0, j = 1 of the operator B(j). Coefficient x±(z) will depends on 
^o^o + r i^i arbitrary constants, where r = dimN(B(0)), r\ = dimN (B (1)) . 
Let us introduce the condition 

(C) Let operator B(j) for j e (0, 1, . . . , iV) has regular points only or singular fredholm points 
ji, . . . ,j u of indices fcj, dimN(Bji) = i = 1, za 

Then by the similar means we can determine the rest of the coefficients #2(2), • • • , xn(z) in 
x(t) from the sequence of difference equations 

n-l 

X n (0, 0)x,(z) + ^(a'(O)) 1 ^ (^(0, 0) - K m (0, 0))^(z + a,) + 
i=i 

+P,(xoW, . . . , x^!^))) = 0, j = 27iV. 

Therefore we have the following 
Lemma 1. 

N 

Let conditions (B) and (C) be fulfilled. Then exists the vector- function x(t) = ^Xj(lnt)f, 

i=0 

such as \\F(x(t))\\E 2 = o(t N ), where operator F is defined with formula (4)- The coefficients 

Xi(\nt) are polynomials of In t of increasing powers and they are smaller than ^2,kj singular 

j 

Fredholm points j G {0, 1, 2, . . . , N} of characteristic operator B (j) . Coefficients Xi(lnt) depends 
on Y^j =0 dim.N(B(j))kj arbitrary constants. 

Remark 2. If B(0) is the Fredholm operator and dimN(B(0)) > 1, then coefficient Xo(lnt) can 
be linear function of lnt and vector-function x(t) will increase unboundly for t — > +0 (briefly, 

x e C mT] . El ))- 
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4. Theorem of existence of continuous parametric solutions 



Since < a-(O) < 1, a^O) = 0, i = l,n - 1, then for any < e < 1 exists T G (0,T] such 

as ma x |a-(t)| < e and sup - 
j=i,n-i,te[0,T'] i=T^=T,te(o,T'] 

Let us introduce the condition 



(D) Let operator K n (t,t) has inverse bounded operator for t G [0, T'\ and iV* is selected to 
have the following inequality fulfilled 



n— 1 

sup e^K^I 
te(o.T') ~ 



K~ {t, t) Ki(t, a t (t)) - K i+1 (t, Oi(t)) 



<q<l. 



C{Ei^E{) 



Lemma 2. 

Let condition (D) be fulfilled. Let in C([ 0i r'];Bi) exists element x(t) such as for t — > +0 

||F(x(t))|k = o(^),iV>J\r. 
TTien equation (3) in C([ 0; t'];Ei) has solution 

x(t)=x(t) + t N *u(t), (22) 
where u(t) is unique and defined with successive approximations. 



Proof. Let us substitute (22) into the equation (4). We get the following integral-functional 
equation for determination of function u(t) 



K n (t, t)u(t) + £ a[(t) N (Ki(t, Oi(t))- 



-K i+1 (t,ai(t)) )u(aii(t))+ 



(*i-i(t) 

Let us introduce the linear operators 



+ E J K ^^ s X t ) u{s)ds + F{x(t))/t N * =0. 



Lu d ^ K-\t, t) £ <(t) (^P) N [m, a t (t)) 



-K i+1 (t,ai(t))ju(ai(t)), 
Ku d ^J2 / K 1 (t,t)K^(t,s)(s/t) N \(s)ds. 

1=1 ,,s 



(23) 



«i-i(t) 
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Then system (23) can be rewritten as follows 

u + (L + K)u = 7 (i), 

where ^(t) = K~ 1 (t,t)F(x(t))/t N * is continuous vector-function. Let us introduce the a Banach 
space X of continuous functions u(t) with ranges in a Banach space E 1 and following norm 

I lull; = max e~ lt \\u(t)\\ El , I > 0. 

II II' < t < T / M v 

Norm of linear functional operator L satisfies the estimate 

||£|U(x-^x) < q < ■ 

because of inequality sup < e < 1 and based on the condition (D) V/ > 0. Moreover for 

t€(0,T'] 

integral operator K the following estimate is fulfilled 

\\K\\ C (x^x) < qi < i - q 

for big enough I. As results, for big enough I > we have 

\\L + K\\c(x->x) < 1, 

i.e. the linear operator L + K is contracting in space X. Therefore the sequence {u n } converges. 
Here u n = -{L + K)u n -i + i(t), u = i(t). 

Theorem 2. (Main Theorem) 

Let conditions (B), (C), (D), /(0) = be fulfilled. Let operator 5(0) has bounded inverse 
operator. Then equation (3) in space C([ 0i t];Ei), < t < T' < T has solution 

x(t) = x(t) +t N *u(t), 

V 

which depends on r iki arbitrary constants. Moreover, element x can be constructed as logarithm- 
ic 

power sum (9), and u(t) is uniquely computed with successive approximations. The asymptotic 
estimate \ \x(t) — x(t)\\ El = 0{t N ") is fulfilled for t — > +0. 

Proof. Let us employ Lemma 1 and take into account conditions of the theorem. Therefore the 
construction of asymptotic approximation x(t) of the desired solution as following logarithm- 
power polynomial 

N 

^Xi(lnt)f 

is possible. And coefficients Xi(\nt) depend on certain number of arbitrary constants. We can 
now apply Lemma 2 and x(t) = x(t) + t N *u{t). Therefore continuous function u(t) can be 
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constructed with successive approximations. Theorem is proved. 



The parametric family of solutions constructed on [0, T"] we can continued with "step 
method" [8, p. 199] on the interval [0,T]. 

If j — is Fredholm point of operator B(j) and dimN(B(0)) > 1 then based on Remark 3 the 
coefficient xo(lni) in the asymptotic x(t) can be linear function of Int. In this case the solution 
x(t) G C(( ,T];Bi) and grow unboundedly as t — > +0. 

Example 3. 

Equation 

t/2 1 t 1 

j J K(y,y 1 )x(s,y 1 )dy 1 ds + J (^J K{y,y l )x(s,y 1 )dy l -2x{s,y) S j ds = g(y)t, 

t/2 

where < t < oo, < y < 1, 1 is eigen value of continuous symmetric kernel K(y, yi) of rank r, 
{0i (y), . . . , 4> r (y)} is corresponding orthonormal system of eigen functions on [0, 1], g(y) G C[ 0i i], 
meets the conditions of Theorem 2. In such case j = is simple singular Fredholm point of 
corresponding characteristic operator B(j). The solution of equation is following 

l 

111 t> ^ f 

x(t, y) = - / Mv^iiy^fiyi) d yi + c i<M2/) + • • • + c r r (y) + x (y), 
i=i ^ 

arbitrary constants, x (y) is the particular solution of Fredholm integral equation 
of the second kind 



1 r 1 

:(y)= / K{y,y l )x{y 1 )dy l - f(y) + J^&(y) / <fyi 

n 1=1 n 



1 

CI?/) = 



An Improvement of the Theorem 2. 

Let {j!,...,^} G N|J{0} are Fredholm points of characteristic operator B(j). Let us 
construct generalized Jordan sets in sense of [19, p. 30] for the characteristic operator B(j). 
Then Theorem 2 can be improved. Indeed let j* is singular Fredholm point of the operator 
B(j). Let elements {(f)f^},i — l,r,l — l,Pi are constructed and satisfy the following equalities: 

= o, 

5(j*)0f ) + J B( i )( J *)(% +i )0; i) = o, 

(24) 

% = l,r, / = I, pi - 1. 
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Let also the following estimate be fulfilled 

i,j=l,r 

where {<pi} r \ i s basis in N(B'(j*)). Then we can follow known theory of Jordan sets of linear 
operators ( readers may here refer to [19, ch.30]) and say that operator B(j) for j = j* has 
complete B(j*) — Jordan set (CJS) (0f^)j=i^j = i^T- We name pi as lengths of Jordan chains 

It is to be noted that CJS exists in j* if j* is singular Fredholm point of index p in sense of 
Definition 3. In this case B^(j*)<f)f^ — 0, % — l,r, I — 0,p — 1, Jordan chains 

are stationary, i.e. they have the same lenght p, Condition (25) become following 

det[< B^4\^>]^ =T7r ^0. 

Also condition (C) satisfied, i.e. we get the result of the Theorem 2. 
Let us relax condition (C): 

(CI) Characteristic operator B(j) for j G (0,1,..., N*) has exactly v singular Fredholm points 
(ji, . . . ,j v ) with complete generalized Jordan sets, the rest of the values of this array are 
regular. 

It is to be noted that in condition (CI) Jordan chains can be nonstationary. 

If condition (CI) is satisfied then an asymptotic approximation x of desired parametric 
family of solution of the equation (4) can be constructed. Indeed let operator B(0) is Fredholm 
operator, {(j)^Y i=l is the basis in N(B(0)), {0^}j = T7 i=Tpl ls corresponding CJS satisfying 
equations (24) and condition (25) for j* = 0. Then the first coefficient Xq(z) of desired 
approximation x, satisfy the difference equation (12) and can be constructed as polynomial 

r Pi 

x (z) = J> £ 4>?^ +1 ~ l + *o, (26) 

i=l 1=1 

where constants and x G E\ to be determined. Substitution (26) in (12) yields 

linear equation wrt x(0) 

r 

B(0)x + CiiB^iO)^ + ^^-^(O)^ + • • • + 5 (1) (0)4 Pl) ) = /'(0). (27) 
i=i 

Because of inequality (25) for j* = we can fine vector (ci,...,c r )' from the condition 
for solvability of the equation (27). Next we can apply Trenogin regularizer ([20], p. 221) 

14 



det 



< B^\f)^ + B^if)^ + ■■■ + B^(f)^\^ 



> 



and construct the solution x of inhomogeneous equation (27) up to the basis . . . , r ) in 
N(B(0)) Similarly, due to the condition (CI) the rest of the coefficients x 2 (z), . . . , xn(z) can be 
computed. Therefore based on the Lemma 2 we got Theorem 2 improved due to the condition 
(C) relaxation: 

Theorem 3. Let conditions (B), (D), (CI), /(0) = be satisfied. Then equation (3) for 
< t < V < T has parametric family of solutions. 

Conclusion 

In case of single equation when E\ = E 2 = R 1 , our method for solution of the difference 
systems solution turn to be known A.O.Gelfond's method (readers may refer to [18], p. 338) 
of construction of particular solutions of inhomogeneous difference equations with polynomial 
right-hand side. We employed the results from functional analysis [15, 19, 20] and the ideas of 
the Gelfond's method have been applied in the theory of the Volterra linear equations of the first 
kind with piece-wise continuous operators acting on a Banach spaces. The constructed method 
can be employed for solution of the class of integral-functional systems with singularities [14, 
15]. If /(0) 7^ then the equation (1) does not have solution in class of continuous functions. 
In this case our method and results from the paper [13] enable constrution of the solutions in 
class of distributions [16, 17, 21]. 
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